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Abstract The approach to study properties of charge-exchange excitations in hot nuclei 
is presented. The approach is based on the extension of the finite rank separable approxi- 
mation for Skyrme interactions to finite temperatures employing the TFD formalism. We 
present the formulae to obtain charge-exchange strength distributions within the Thermal 
q ■ Quasiparticle Random Phase Approximation (TQRPA). 
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The properties of nuclei at finite temperatures appear to be interesting for many reasons. 
Currently, one of the most popular fields to apply the nuclear theory at finite temperatures 
is the astrophysical studies. In the astrophysical context the theory is used to calculate 
thermal modifications of spin-isospin transition distributions over nuclear spectra, as they 
play an important role in weak- interaction mediated reactions in stellar environment. 

At first, temperature effects were introduced into the calculations via quite straight- 
forward way taking into account a possibility of thermal feeding of nuclear excited states 
which energies were calculated via independent-particle shell model or taken from exper- 
imental compilations [I]. Later this approach was developed and refined employing the 
large scale shell- model calculations of low-lying excited states j2]. Within this approach 
a coupling with high- lying resonant states like, e.g., the Gamow- Teller resonance, was 
\Q • considered in the framework of the Brink hypothesis. 

The other approach to calculate the rates of weak-interaction mediated processes 
at finite temperatures employs the thermal quasiparticle random phase approximation 
(TQRPA). For the first time it was used in [3] to study the electron captures on neutron- 
rich nuclei. Recently, there appeared several papers |H |SJ E] where the authors applied 
TQRPA (or TRPA) with different nuclear Hamiltonians. 

In Ref. [5], the TQRPA was applied in a framework of a general context of thermo- 
field dynamics (TFD) [TJ |8] and with the Hamiltonian consisting of schematic separable 
effective interactions. It is well known that schematic separable interactions allow one to 
avoid many of computational difficulties and limitations which are inherent for "realistic" 
residual interactions. However, the cost is a lack of selfconsistency and limited predictive 
power of the calculations for nuclei far from the valley of stability. Thus, it seems reason- 
able to use in calculations for hot nuclei a finite rank separable approximation elaborated 
for the residual forces based on an effective interaction of Skyrme type [T0~| [TT| 112] . This 
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enables to combine the advantages of consistency (the mean field and the residual inter- 
action of RPA are determined from the same effective interaction) with a computational 
simplicity (the size of the RPA problem does not increase with increasing configuration 
space) . 

Here we present the formulae which extend the approach of Refs. [10] to finite tempera- 
tures. The TFD formalism presented in [HI E] and adapted to nuclear structure problems 
in [13] is used. Having in mind further application of the formalism to charge-exchange 
excitations we consider only the isovector part of the residual particle-hole interaction. 



2 Separabelization of residual interactions 

The exact p-h residual interaction V res corresponding to the Skyrme force can be obtained 
as the second derivative of the energy density functional with respect to the particle den- 
sity. In Ref. JTU], V res was approximated by its Landau-Migdal form. For the Skyrme 
interactions all the Landau parameters with I > 1 are zero. Moreover, the Landau pa- 
rameters with the / = 1 terms are neglected. Therefore, the residual interaction in the 
isovector channel has the following form: 

V res (ri, r 2 ) = N^[F>( ri ) + G^O^V^'V 2 ^ - r 2 ). (1) 

where and 7~W are the spin and isospin operators, and iVo = 2kpm* /n 2 h 2 with kp and 
m* standing for the Fermi momentum and nucleon effective mass in nuclear matter. Fq, 
G' are functions of the coordinate r. Their expressions in terms of the Skyrme parameters 
can be found in Ref. [14j . 

Following the method presented in [101 [11] the residual interaction is reduced to a 
finite rank separable form 

V res = V M + V SM , (2) 



where 
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VSM — -2 ^ ^ ^2 K G } '■ ^LJM^LJM : > ( 3 ) 
JM L=J;J±l n=l 



N^^'l (4) 



and 

Here R is a large enough cutoff radius for a iV-point integration Gauss formula with 
abscissas r k and weights Uk [ID] . 

The operators entering the normal product in Eq. ([3]) are defined as follows: 



Mffl = {-l) J - M J- 1 Y, ti" m ^ ~ m P \ JM )ff2i- l ) jV '^ a U n ^ 

Sfl} = {-l) J ~ M J' 1 £ (^m^-m p |JM)^(-l)^-^at nmn% ^, (5) 
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where fj^j is the single-particle matrix elements of the multipole operators, 

fiS = u jn (r k )u jp (r k ){j n \\i J Yj\\j p ), (6) 
and where gfj^ is the single-particle matrix elements of the spin-multipole operators, 

9jtj p k) = u u ( r k)u jp (r k ) (j n || i J Tj || j p ) . (7) 

In the above equations, (j n ||i J Yj||j p ) is the reduced matrix element of the spherical har- 
monic Yjm, J = v2J+l, Tljm = \Xj x a] J M . The radial wave functions %(r) are related 
to the HF single-particle wave functions 

0i,m(l) = ^3% ( (fl,a 1 ). (8) 

An interaction in the particle-particle channel has the surface peaked density- dependent 
zero-range shape 

v pair (n, r 2 ) = v (i - v (r^r) ) 5 ^ - r 2)- ( 9 ) 

Here p(r) is the particle density in coordinate space, p is equal to the nuclear saturation 
density, Vo, rj and a are the parameters fixed to reproduce the odd-even mass difference of 
nuclei in the study region. The interaction (Q is responsible for the pairing correlations 
as well. 



3 TBCS and TQRPA equations at T ^ 
3.1 Fundamentals of the thermo-field dynamics 

Thermo-field dynamics [3 El [9] is a real-time formalism for treating thermal effects in 
quantum field theory and many-body theories. In TFD, the thermal average of a given 
operator A is calculated as the expectation value in a specially constructed, temperature- 
dependent state |0(T)) which is termed the thermal vacuum. This expectation value 
is equal to the usual grand canonical average of A. In this sense, the thermal vacuum 
describes the thermal equilibrium of the system. 

To construct the state |0(T)), a formal doubling of the system degrees of freedom is 
introduced. In TFD, a tilde conjugate operator A - acting in the independent Hilbert 
space - is associated with A, in accordance with properly formulated tilde conjugation 
rules |Zl El E] 

(A 1 A 2 y= MM, {cxAx + c 2 A 2 r= c\M + c* 2 A 2 , 

(A^=(A)\ (AT=A, (10) 

where A± and A 2 stand for any operators and c% and c 2 are c-numbers. The asterisk 
denotes the complex conjugate. It is assumed that any bosonic (fermionic) operator with 
tilde commutes (anticommutes) with all bosonic (fermionic) operators without tilde. 
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Let H be the Hamiltonian of the system. In the doubled Hilbert space the thermal 
vacuum is defined as the zero-energy eigenstate of the so-called thermal Hamiltonian 
H = H — H and it satisfies the thermal state condition [7J [SJ |H] 

A\0(T)) = ae n / 2T A^\0(T)} } (11) 

where a = 1 for bosonic A and a = i for fermionic A. The thermal state condition is 
one of the most fundamental relations in TFD. The celebrated Kubo-Martin- Schwinger 
condition, which is one of the basic axioms in the c* -algebra formalism for statistical 
mechanics, is a result of the thermal state condition in TFD. Furthermore, from (JTTJ) we 
can see immediately that, in TFD, there always exists a certain combination of A and A^ 
which annihilates the thermal vacuum. That mixing is promoted by a specific canonical 
transformation called the thermal Bogoliubov transformation [S]. The temperature 
dependence comes from the transformation parameters. 

The important point is that in the doubled Hilbert space the time-translation operator 
is the thermal Hamiltonian % = H — H . This means that the excitations of the thermal 
system are obtained by the diagonalization of "H. The existence of the thermal vacuum 
annihilation operators allows for straightforward extensions of different zero-temperature 
approximations to diagonalize the thermal Hamiltonian. As follows from the definition of 
% each of its eigenstates with positive energy has the counterpart - the tilde-conjugate 
eigenstate - with negative but the same absolute energy value. Transitions from the ther- 
mal vacuum to positive (non-tilde) energy states correspond to excitation of the system, 
while transitions to negative (tilde) energy states - to deexcitation. 

3.2 Equations for pairing correlations at T^O 

Applying the TFD formalism we suppose that the nuclear proton and neutron mean 
fields are already produced according the Hartree-Fock procedure and our Hamiltonian 
H consists of the mean fields and separabelized residual interactions presented in the 
previous Section 2. In particular, it means that we ignore the influence of temperature 
on the nuclear mean field. 

At first, following the TFD prescription, we double the original nuclear degrees of 
freedom introducing the so-called tilde creation and annihilation operators aj m ,a., m and 
construct the thermal Hamiltonian 

U = H-H 

Then we transform % to the thermal quasiparticle representation by means of the two 
canonical transformations. The first one is the standard Bogoliubov transformation to 
quasiparticle operators 

a jm ~ U j a jm V j a Jm 

Ojm = U Am + K 2 + V ) = !)> ( 12 ) 

where a-j^ = (— l)- 7 ' -m ay_ m . The second transformation is the thermal Bogoliubov trans- 
formation IB]. It mixes the quasiparticle and tilde quasiparticle operators, thus pro- 
ducing thermal quasiparticle operators and their tilde partners (3j m , /3 jm , /?j m , /?- m 

fijm = X j a ]m ~ Wj a jm j 

Pjm = X i^) m + Wj(Xjm , {X* + Vj = 1) • (13) 
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The coefficients of both the transformations are determined by the diagonalization 
of the sum of single-particle and pairing parts of H and additional demand for the BCS 
thermal vacuum to obey the thermal state condition [T5] . At the end, we get the following 
equations for the coefficient Uj,Vj and Xj,yj\ 

«i = ^f 1 -— y^Hi-^) 172 , (14) 



V2 

Vi = 



exp 



Ej- 


A 















xj=(l-y*) 1/2 , (15) 



where Ej = y (Ej — A) 2 + A 2 - is a quasiparticle energy. The pairing gaps Aj and the 
chemical potential A are the solutions to the finite-temperature BCS equations 



j' » 

iV = ^(2j + l)(^ 2 + n>|), (16) 

3 

where N is the number of neutrons or protons in a nucleus, and two-body matrix element 
Go(jjj'j') is given by the J = particle-particle matrix element of the interaction V pa i r 

Go(jJj'j') = (jj\V pair \j'j') 00 . (17) 

Now, the sum of the single-particle and the pairing parts of the Hamiltonian 7-L becomes 
diagonal 

^BCS ^ ^ ^ ^ (j^jmPjm PjvaPjm^ ■ 
r=p,n jm 

Here implies a summation over proton or neutron single-particle states only. The 
Hamiltonian "Hbcs describes a system of noninteracting non-tilde and tilde thermal quasi- 
particles with energies Ej and — Ej, respectively. The vacuum of thermal quasiparticles, 
|0(T); qp), is the thermal vacuum in BCS approximation. 



3.3 Thermal QRPA equations 

The transformations (|12p and (fT51) are applied to the whole thermal nuclear Hamiltonian 
H. At the next step we approximately diagonalize H within the Thermal Quasiparticle 
Random Phase Approximation 

H ~ HtQRPA = ^2 Uji ^JMiQjMi — Q^jMiQjMi)- (18) 
JMi 

Here the thermal (charge-exchange) phonon creation operator is defined as a linear su- 
perposition of the proton-neutron thermal quasiparticle pair creation and annihilation 
operators 

Jnjp 

+ €^M J M + ^IPMi + i^Mi + ^XlPMi), (is) 
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and [ ] J M denotes the coupling of single-particle angular momenta j n , j p to total angular mo- 
mentum J. The tilde conjugate thermal phonon operator Q\ Mi can be obtained from ffT9]) 
by applying tilde conjugation rules f lTO]) . Now the thermal equilibrium state is treated as 
the vacuum |0(T);ph) for thermal phonon annihilation operators and it obeys the ther- 
mal state condition (TTTT) . The excited thermal one-phonon states are Q^ JMi \0(T); ph) and 

The thermal phonon operators are assumed to commute as bosonic operators, that is, 
[Q jmvQj 1 m'A = &JJ'$MM'$ii> ■ This assumption imposes the normalization condition on 
the phonon amplitudes 

V(V ; ' il>F- + ■ ti Ji '■ ■ if' 1 n' 1 ' +n Ji n Ji '- 

/ / \' Jnjp ' Jnjp 'Jnjp' Jnjp 'Jnjp 'Jnjp 'Jnjp 'Jnjp 

Jnjp 

^jnjp^jpjn ^jpjn^jnjp ^jnjp^jnjp ^jnjp^jnjp^ ( ^) 

Furthermore, additional constraints on the amplitudes come from the thermal state con- 
dition. Namely, putting A = [a^ajj^ in (fTTT) we get (uji > 0) 

(VjpTPjljp + yjnVjMdp) = eXP (lf ) (V^^ + yjnVjp^JnJp)^ 

(Vip^jp + yjnVjp(Pjt jp ) = exp (-^) (Vj P <f>% jp + UjJljSfj;)- (21) 
and for A = [a] n a^]jv/ we have 

^jnVjp^jnjp + Vjn X jpij n j p ) = 6X P (^J^) ( X JnV jp^ j„jp + V jn X jpV j n j p ) , 

( X jnyjpVj„j p + yjn X jpCj„j p ) — eX P ( _ ^JnVjp^jnjp ^ & 'jn X jpVj„j p ) • (22) 

To find the energy and the structure of thermal phonons we apply the variational 
principle, i.e. we minimize the expectation value of "H over the thermal one-phonon 
state under the constraints fl20l) . As a result we get the TQRPA eigenvalue equations 
for the amplitudes and the energy eigenvalues uji. In contrast to the zero-temperature 
case, the negative eigenvalues of TQRPA matrix have physical meaning and they can be 
interpreted as the excitation energies of tilde thermal one-phonon states relative to the 
thermal vacuum. Besides, each eigenvalue is twice degenerate so that %tqrpa ([18]) is 
invariant under the thermal Bogoliubov transformation 

QjMi ~~ XjiQ JMi — YjiQjMi, Q\ Mi — > XjiQ\ Mi — YjiQjMi (23) 

with Xj { — Yj { = 1. To find the structure of thermal phonons unambiguously we demand 
that the constraints (12 1 |I22[) are valid. Only in this case the vacuum of thermal phonons 
is the thermal vacuum in the TQRPA approximation. 

Using the separable form of the residual interaction one can reduce remarkably the 
dimensions of the corresponding TQRPA matrixes. To do this we introduce a vector 
by its components: 

D ± n = J2€t{ u %p[^^^ 

Jnjp 

~^ V Jnjp [ X JnV jp(Vj n j p =t ij n j p ) =t yjn X jp(Vj„j p =^ t,j n j p )] j*) (24) 
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where uf 1 } — Uj Vj ±i>,- Uj and v^} = Uj Uj ±Vj Vj . Phonon amplitudes are functions 

Jnjp Jn Jp Jn Jp Jnjp Jn Jp Jn Jp r 

of the vectors 

2N 



w 




7-2 
J 


<t>, 

T / 


L- 

Jnjp 


e {+) tout- 

b JnJp + U ^ 


/ \ 

V 


Ji 


j — o 
J Z 
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)..= 

Jnjp 


£ jnjp ± U -Ji 




Ji 


J' 2 


V 


= 

Jnjp 






Ji 




u 


Jnjp 


4i ± ^ 



- (vj p Xm - yjnVjpYxi) E 4£" (B) ± ^U) . 

n=l 
2 AT 

n=l 

- v*^0 E 42^ ± ^ 



Jnjp 



n=l 
2N 



Jin v (+) 

Jnjp J ' 



- (.o,//^A,--//,, a .v A; ) e 42« (n) (^>U =f ^-Ht) , (25) 



n=l 



where 



-1/2 



L exp(^)-i]'-'", x Jl = (i + r A ^/ 2 . 

So, the TQRPA equations are reduced to the set of equations for D± n : 

Mi-y m 2 



m 2 M 3 -y 



D. 



0. 



(26) 



(27) 



The matrix elements of the 2N x 2N matrices j\4p have the following expressions 



(ri) f c ( + ) C.W ^2 

K ,(Jn)JJn') J £ j n jp\ U jnjp> 



MT > = l^rfflyj (i-yl-yl) 

J-w 72 / / Jnjp Jnjp 1 / (+) \2 2 J™ jp' 



eH ( v ^ ] ) 2 

C -JnJp\ U JnJp> / 2 _ 2 

/--.(-) ^2 _ . ,2 g 3p> ( 



MT 



K 



(n') 



J 2 



--E424 



(Jn)i(Jn') J a jnjp a jnj P 



Jnjp 



* I ^(+) ) 

"Jnjp' 



— UJ 



i ~ yjn ~ Vjp) ~ , (-) N9 ^(y^n — y? P ) 



(4i) 2 -^ 2 



Ji 



(28) 



where e^j = £j n ± Ej p . Its solution requires to compute the determinant 

V A^ 2 M 3 -\I ) 
and we find the eigenvalues of the TQRPA equations. 

4 Charge-exchange transition probabilities 

Charge-exchange transition probabilities (transition strengths) from the thermal vacuum 
to thermal one-phonon states are given by the squared reduced matrix elements of the 
corresponding transition operator: 



*L ±} = 



<Q JM ,||£jMl|0(T);ph> 

<Qjtfill42ll°C0;ph> 



(29) 
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Hereinafter the symbol (— ) labels the n — > p transition operators, and the symbol (+) 
labels the p — )■ n transition operators. The explicit expressions for $ g° and fc^f are the 
following: 

Jnjp 
Jnjp 

= (E(-i) 5Wn+J 4 -) (ipin)^(inip;Ji)) 2 , (30) 

where c^^OpM-Mp)) is a reduced single-particle matrix element of the transition operator 

d { \\ 3 p{n) 3 n(p)) = (jn(p)\\D^ \\j p (n)), (31) 

functions Qi^ijnjp] Ji) are linear combinations of the phonon amplitudes (|25|) : 

Jt ~ U jn U jp( X j n yjpVj n j p +yjn X jp£j n j p ) ~ V jn V jp(yj n X jpVj n j p + X in%'pC 7 „ Jp ) ) 

fi 2 (jnj P ; Ai) = v jn u jp (x jn x jp ^ jp +y jn y jp ^ jp ) +u jn v 3p (y jn y jp ^ jp +x jn x jp (j)^ jp ) - 

+ v j^ P {x jn y jp 7]^ Jp +y jn x^ jp ) +u jn u jp (y jn x jp ?}^ jp +x jn y jp ^ jp ), (32) 

and f2 1>2 (jnjp; ^) result from fli^ijnjp] Ji) by changing non-tilde phonon amplitudes by 
their tilde partners and vise versa. 

The excitation energies with respect to the thermal equilibrium state of the parent 
nucleus are 

£T = Uj . zp AA (33) 

for non-tilde phonon states, and 

E T = -tojt T AA (34) 

for tilde phonon states. Here AA = A n — A p is the difference between the neutron and the 
proton chemical potentials. 

Expressions (130]) and (1331 134"|) determine charge exchange strength distribution in a hot 
nucleus within the TQRPA. Note that at finite temperatures some amount of transition 
strength is always located in the region of negative transition energies. 
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